Spectral Functions of One-dimensional Models of Correlated Electrons 
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Using the Ogata-Shiba wave function, the spectral functions of the one- dimensional (ID) infinite 
U Hubbard model are calculated for various concentrations. It is shown that the "shadow band" 
feature due to 2/cf fluctuations becomes more intense close to half-filling. Comparing these results 
with exact diagonalization data obtained on finite clusters for the finite U Hubbard model and for 
the t — J model, it is also shown that this feature remains well-defined for physically reasonable 
values of the parameters (U/t ~ 10, J/t — 0.4). The "shadow" structure in the spectral functions 
should thus be observable in angle-resolved photoemission experiments for a variety of quasi-one 
dimensional compounds. 
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Spectral functions are very useful to understand the 
electronic structure of solids. They are defined by 

A(k,u>) = £ \(f,N+l\clJO,N)\ 2 8(uj-E? , 

B(k,u;)=J2\(f,N-l\c k jO,N)\ 2 S(u J -E^ + Ef~ 1 ), 
f," 

in the standard notation, and in principle they can be 
measured in angular-resolved inverse photoemission and 
photoemission experiments, respectively. The effect of 
electron-electron interactions on the spectral functions 
of the three-dimensional Coulomb gas has been investi- 
gated in much detail several years ago Q . With respect 
to the simple S(u> — £&) structure of the non- interacting 
case, the main differences are: i) A shift of the energy 
of the quasi-particle band; ii) A broadening of the quasi- 
particle peak; iii) A renormalization of the weight of the 
quasi-particle peak, compensated by the appearance of 
an incoherent background; iv) The presence of another 
band - a plasmon band - due to the long-range nature of 
the Coulomb potential. For a short-range repulsion, such 
as the on-site interaction of the Hubbard model, there is 
no plasmon band, and the spectral function is expected 
to have only one well-defined feature, the quasi-particle 
band, on top of an incoherent background. The validity 
of this simple picture in lower dimensional systems is cur- 
rently under intense discussion. In two dimensions(2D), 
there are at least two important issues: 1. the actual ex- 
istence of a quasiparticle peak, and 2. the generation by 
antiferromagnetic short-range fluctuations close to half- 
filling of additional well-defined features in the spectral 
functions, called generically "shadow bands" The ab- 
sence of exact results in 2D for correlated models makes 
the interpretation of the experimental results quite diffi- 
cult. For high-T c cuprates, these issues are still contro- 
versial 101. 



The situation is much more favorable in ID. First of 
all, it is known that there is no quasi-particle peak in 
the spectral function but two divergences due to spin- 
charge separation Q. Besides, using the Ogata-Shiba 
wave- function exact results have been recently ob- 
tained H for the spectral functions of the infinite U limit 
of the Hubbard model defined by: 



or, equivalently, for the J 
defined by: 



U^n lA n lA , (1) 



limit of the t — J model 



i,fT 



(2) 



where c are the usual hole operators, and the rest of the 
notation is standard. Combining these results with nu- 
merical simulations should allow us to reach a precise 
picture of the spectral functions of one-dimensional sys- 
tems. 

In this paper we will concentrate our efforts on three 
aspects of the spectral functions of ID systems: 1. The 
first one is the extra weight created near the Fermi en- 
ergy close to half-filling by short-range antiferromagnetic 
correlations. It is this feature that is usually referred 
to as "shadow band" in the literature. 2. The second 
issue presented here is that the "shadow" weight is ac- 
tually part of an extended band that can be followed in 
the whole Brillouin zone. 3. The third point is that this 
extended band still exists away from half-filling due to 
diverging 2k p fluctuations, although its weight near the 
Fermi level diminishes rapidly with hole doping. In por- 
tions of the paper we will refer to all these features with 
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the common language of "shadow band" , but it must be 
clear to the reader that in doing so we are extending the 
usual meaning of this term to include all the extra inter- 
esting features caused by 2Uf fluctuations, both close and 
away from the Fermi energy, that appear in the spectral 
functions of ID systems. 

The presence of shadow bands in ID systems was first 
revealed in studies of the spectral function correspond- 
ing to a hole injected into one dimensional Heisenberg 
models at half- filling Q|. They were identified as poten- 
tially observable features in photoemission experiments 
for CuGeC>3 and other undoped compounds where anti- 
ferromagnetic fluctuations are important. Other related 
aspects were first discussed in Ref. || where a calculation 
of the spectral functions in the infinite U limit showed un- 
ambiguously that, in addition to the charge and spin fea- 
tures typical of a Luttinger liquid, there is a well-defined 
structure whose dispersion can be followed in the whole 
Brillouin zone. This accumulation of weight in the spec- 
tral functions originates from diverging spin fluctuations 
H at 2kp, the ID equivalent of the diverging fluctua- 
tions at (Tv/a,n/a) in the 2D case at half- filling, and by 
analogy it was also called a shadow band. 

A clear identification of these features in photoemission 
experimental results for quasi-one dimensional systems 
is still lacking. On the theoretical side, progress can be 
made along two lines: First, by estimating the evolution 
of the shadow band features, in particular their intensity, 
with electron concentration to determine which systems 
are a priori the best candidates to observe the effect. 
Second, by extending the calculations to more realistic 
parameters, i.e. U < +oc or J > 0, to investigate if the 
band remains a well-defined structure when the on-site 
repulsion is no longer infinite. 

To address the first point, we have calculated the spec- 
tral functions for an infinite repulsion along the lines of 
Ref. H for various concentrations (9). Some results for 
B(k,uj) are shown in Figure 1. The evolution between 
n = 1/2 and n = 1 is smooth. When 3kp is larger than 
Tt/a, i.e. for n > 2/3, the band at finite binding energy 
discussed in Ref. |6| at n = 1/2 is now folded back when 
it hits the zone boundary and it crosses the Fermi level at 
2n / a — Zkp. This detail is very important since now non- 
trivial structure in the spectral functions appears near 
the Fermi energy. The rest of the shadow band evolves 
smoothly with density becoming more intense close to 
half-filling. Note, however, that the effect is likely less 
dramatic as expected in 2D systems for realistic values 
of U/t. This is due to the fact that in ID there is always 
perfect nesting, independent of the band-filling, since the 
Fermi surface consists of only two points. The ID best 
candidates to show these shadow band effects are com- 
pounds with a concentration close to 1, as anticipated in 
Ref. (7). However, quarter-filled systems should not be 
disregarded since the portion of the band that appears at 
energies between — 2t and —t for < k < ir/2a could be 
experimentally detectable, unless the usual large spuri- 
ous backgrounds that contaminate photoemission results 



as we move away from the Fermi level hide the signal in 
the noise. 




FIG. 1. Spectral function B(k,u) of the U — > +00 Hub- 
bard model calculated with the help of the Ogata-Shiba 
wave-function, a) n = 3/4 (kF = 3ir/8), 200 sites; b) n = 7/8 
(k F = 7tt/16), 176 sites; c) n = 1 (kp = tt/2), 190 sites. 

The second point, i.e. calculations for realistic pa- 
rameters, is more difficult to address. Even for models 
that are integrable, such as the Hubbard model or the 
super-symmetric t — J model (J/t = 2), a calculation 
of the spectral functions for large systems remains an 
open problem. The only available methods are numer- 
ical |u^]. Relatively large systems can be handled with 
Monte Carlo simulations |0|, but the results require an 
analytic continuation which is difficult to control. Thus 
exact diagonalizations of finite clusters seem to be more 
appropriate. Such calculations have already been pre- 
sented for ID systems |7|, |l2| , p^| . The main limitation of 
these simulations is the size of the systems - typically 16 
to 20 sites close to half-filling. The spectral functions ob- 
tained with the exact diagonalization method consist of 
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a collection of discrete 5 functions, and they are difficult 
to interpret without an underlying theory. 

The originality of the present work is to use the es- 
sentially exact results obtained for infinite U with the 
Ogata-Shiba wave function to interpret the results for fi- 
nite clusters. To illustrate how this works, we have com- 
pared in Fig. 2 the "exact" results for n = 7/8 (Fig. 2a) 
with the results obtained for the same model on a 16 site 
cluster using exact diagonalizations (Fig. 2b). Details 
about the numerical method can be found in Ref. [^0[ . 
Clearly, the calculation on 16 sites is not able to repro- 
duce a continuum: The true spectral function has a con- 
tinuum between the peaks corresponding to the charge 
and spin excitations of the Luttinger liquid Q], whereas 
the spectral function of Fig. 2b has several peaks in the 
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tinger liquid features is not possible on the basis of ex- 
act diagonalizations of small clusters, unless a variety of 
boundary conditions is used and a finite size scaling anal- 
ysis is carried out. In spite of this problem, the shadow 
band is clearly visible on the data of Fig. 2b, and both its 
intensity and its dispersion are in good agreement with 
the exact results. The continuum that exists between the 
Van Hove singularity at — It and the shadow band around 
k = it is again poorly described on the finite cluster, but 
the main band and the shadow band can be clearly iden- 
tified. Thus we believe that discussing the intensity of 
the shadow band on the basis of exact diagonalization 
results is meaningful. 
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FIG. 2. Spectral functions of the U — > +00 Hubbard model 
at n = 7/8 for some selected momenta. Solid lines: B(k,u>); 
Dashed lines: A(k,u). a) "Exact" results obtained with the 
Ogata-Shiba wave-function for 176 sites. The rapid oscilla- 
tions are finite-size effects; b) Numerical results obtained on 
a 16 site cluster using exact diagonalizations and periodic 
boundary conditions. The 8 peaks have been given an arbi- 
trary broadening of 0.03t. 

same energy interval. This is best observed for k = 0, 
where the continuum is replaced by a collection of 4 
peaks. This proves that detecting the position of the 
spin and charge excitations corresponding to the Lut- 



FIG. 3. Spectral functions for the Hubbard and t — J model 
at n = 3/4 obtained on a 16 site cluster using exact di- 
agonalizations and antiperiodic boundary conditions. Solid 
lines: B(k, u); Dashed lines: A(k, w). a) Hubbard model with 
U/t = 10; b) t - J model with J/t = 0.4. 

Let us now turn to more realistic parameters. To make 
our analysis below more convincing, we have calculated 
the spectral functions for both the Hubbard model and 
the t — J model, and we have compared the results us- 
ing the relation J = At 2 /U . Some results obtained at 
n = 3/4 for U/t = 10 and J/t = 0.4 are depicted on Fig. 
3a and 3b, respectively. The results are roughly the same. 
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In agreement with what was stated before, we will not 
try to interpret the details of the fine structure around 
w/t = for small momentum or around u/t = —2 for 
large momentum. What is clearly observable, although 
it carries a small weight, are "shadow" features located at 
the same position as in the infinite U case. Its intensity is 
larger in the t— J case than in the Hubbard case, and this 
difference grows as U becomes smaller (since the agree- 
ment between Hubbard and t — J models is expected only 
at large U/t). Increasing J/t does not suppress the spec- 
tral weight around the shadow band but simply makes 
the feature broader, while decreasing U/t suppresses the 
intensity of the shadow band quite dramatically close to 
the zone boundary. The spectral weight which is lost ap- 
pears around u> — U, a region not shown on Fig. 3b for 
clarity. This likely arises from corrections of order 1/U 
that are not included in the t—J model |Q. However, we 
have checked that, for parameters as small as U/t = 5, 
the shadow band remains quite intense for small wave 
vectors. 
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FIG. 4. Spectral functions for the t—J model at n — 7/8 
obtained on a 16 site cluster using exact diagonalizations and 
periodic boundary conditions. Solid lines: B(k,u); Dashed 
lines: A(k,ui). 

Upon approaching half-filling, the shadow band be- 
comes quite intense. For instance, we show in Fig. 4 the 
spectral functions for J/t = 0.4 at n = 7/8. A photoc- 
mission experiment performed on a compound described 
by such parameters should be able to detect the portion 
of the shadow band the closest to the Fermi energy, and 
actually follow it in a large part of the Brillouin zone. 

In conclusion, we have presented a detailed study of the 
"shadow" band that appears in the spectral functions of 
ID systems due to diverging 2fcp fluctuations B. First, 
we have studied its evolution with electron concentration. 
As anticipated in Ref . , its weight near the Fermi level 
is considerably enhanced close to half-filling. Then, we 
have shown that this feature is quite robust for realistic 
values of U/t (or increasing J/t ). This is true even away 
from half-filling specially regarding the weight located at 
a finite binding energy 0. We believe that it should be 



possible to detect this band experimentally either close 
to the Fermi level in nearly half-filled systems or at ener- 
gies of the order of the band width for systems far from 
half-filling, and it is our hope that the present study will 
encourage experimentalists to search for this novel fea- 
ture of the spectral functions of ID correlated electron 
systems. 
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